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Introduction. The use of probabilistic analysis is important when the input data are random, that leads to stochastic 
results. This paper describes the integration of a probabilistic design strategy of the solid and hollow stems implanted in 
a proximal femur in order to compare their advantages. The used hollow stem is called “Improved Austin-Moore” 
(IAM) model. 

Materials and Methods. Probabilistic methods allow variations in factors which control the biomechanical effects of the 
implanted femur to be taken into account while determining its performance. Different material properties were 
generated randomly using Monte Carlo simulation (MCS). Monte Carlo sampling techniques were applied, and 
different von Mises stresses of the layers (bone and metal) were chosen as a performance indicator. 

Results. A simple 2D implant-bone study of solid and IAM stem design was carried out with a high level of confidence, 
99.87%, which corresponds to a target reliability index with regard to statistical uncertainties. The probabilistic design 
results show that the input and output parameters for the IAM stem are highly correlated relative to those for the solid 
stem. 

Discussion and Conclusions. The sensitivity analysis shows that the input parameters for the IAM stem play a much 
larger part in the output parameters relative to the solid stem. The IAM stem is much more advantageous than the solid 
stem which causes an increase in the performance of the hip prosthesis. 
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1. Introduction 


In probabilistic studies on Total Hip Replacement (THR) models, several important random parameters may 
affect the hip performance, such as the geometrical description of the bone, the implant and the cement, the material 
properties of different layers, the magnitude and direction of the applied loads, and the position of the bone-implant 
[1-2]. In order to improve the stem design, several points have been considered. At the first stage, the rectangular 
section stem 1s considered more reliable than the circular section stem [3]. The rectangular section stem then provides a 
secure diaphyseal press-fit in the frontal plane of the femoral canal. This way it allows excellent rotational stability and 
improves the primary mechanical fixation [1]. In addition, wear as mentioned in Kharmanda and Albashi [4] can be 
reduced. At the second stage, the shouldered stem is considered more reliable than the non-shouldered one. It has been 
shown in the numerical results of Kharmanda, et al. [5], that, when considering the shouldered stem, there is an 
excellent decrease of the maximum values on von Mises stress in different bone layers compared to the non-shouldered 
one. At the third stage, the hollow stems are considered. In Kharmanda [6], a multi-objective structural optimization 
strategy has been integrated into Austin-Moore prosthesis in order to improve its performance. The resulting model was 
called Improved Austin-Moore (IAM). The topology optimization was considered as a conceptual design stage to 
sketch the IAM stem according to the daily loading cases. The shape optimization presented the detailed design stage 
considering several objectives. A new multiplicative formulation has been proposed as a performance scale to define the 
best compromise between several requirements. In this paper, a probabilistic design strategy is applied on two kinds of 
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stems (solid and IAM stems) to determine the mechanical effects, the response sensitivities with respect to input 
parameters, and the correlation of the used material properties with different output parameters. A numerical application 
on a 2D problem is carried out to show the advantages of the IAM stem relative to the solid one. 

2. Material and Methods 

2.1 Model description and material properties 

Fig. la and b show 2D sections for the studied solid and IAM stems implanted in bone layers (cortical bone 
and cancellous bone). The bone materials are generally anisotropic [7]. In the current study, the material properties of 
bone have been for simplicity considered as linearly elastic and isotropic. In Fig. 1, the cortical bone material is 
considered to be homogeneous and isotropic with Young’s modulus £ =17 GPa and Poisson’s ratio v = 0.33. The 
corresponding number of elements for the cortical region is 605 elements for the solid stem model, while it is 545 
elements for the IAM stem model. The cancellous bone material is also considered to be homogeneous and isotropic 
with Young’s modulus EF =386 MPa and Poisson’s ratio v = 0.33 [8]. The corresponding number of elements for the 
cancellous region is 417 elements for the solid stem model, while it is 418 elements for the IAM stem model. The 
modulus of titanium alloy of stem is considered to be: L=110GPa with Poisson’s ratio: v=0.3 [9]. The 
corresponding number of elements for the metal region is 529 elements for the solid stem model, while it is 861 
elements for the [AM stem model. The used element is PLANE82 (8-node, nonlinear). The number of the total nodes is 
5048 nodes for the solid stem model, while it is 6094 nodes for the [AM stem model. 
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Fig 1. Geometry models for: (a) solid stem and; (6) IAM one 


2.2 Boundary conditions 
The loading cases are considered as daily loading conditions [10]: one-legged stance (L1), extreme ranges of 
motion of abduction (L2), and adduction (L3) as shown in Fig. 2. 
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Fig. 2. Loading cases: a) One-legged stance case (L1); 5) Abduction case (L2); and c) Adduction case (L3) 


Table 1 presents different components of the three loading cases (L1, L2 and L3) which are illustrated in Fig. 1. 
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Table 1 


Force components 


Loading Body Forces 
Case 


F? [N] F? [N] F” [N] F” [N] 


x y Jy 


It 942.41 -2116.68 -330.04 620.71 
2 -299.7] -1118.54 48.85 347.58 
3 1283.35 -865.63 -268.43 383.36 


According to the results of Kharmanda [6], the third loading case (L3) leads to the highest von Mises stress 
values in the bone tissues. Thus, in this work, a probabilistic design strategy is performed considering the third loading 


Xx 





case L3. 
2.3 Probabilistic design strategy 


2.3.1 Failure probability 

The design of structures and the prediction of their good performance lead to the confirmation of a certain 
number of rules depending on the knowledge of physical and mechanical experience of designers. Thus, there are 
several events leading to a failure mode. It is the objective to evaluate the probability of failure corresponding to the 
occurrence of critical failure scenarios. The probability of failure is then given by: 


P,=Pr[G(xy)]= [| f(ydy,...dy, (1) 


G(x,y)<0 
where G(x, y) 1s a limit state surface (curve) that is defined by the condition of good structural performance [11]. The 
limit state surface (curve) G(x, y)=0 is located between the failure part G(x, y) <0 and the safety part G(x,y)>0. 
F,(y) 1s the density function of the random parameter Y . Several techniques can be used to compute the probability of 


failure. Monte Carlo Simulation (MCS) is the most conservative probability technique [12]. Here, the parameters are 
sampled at random from their underlying distributions and the probability of failure which is estimated by solving the 
model repeatedly. For all time, the MCS gives the correct solution if a sufficient number of trials is computed, but it 
requires a high computing time. In order to decrease this expense, and avoid overlapping of samples and/or lack of 
samples in some regions of the domain, another technique called Latin Hypercube Sampling (LHS) can be applied [13]. 
The MCS based approach divides the sample space into separate intervals with equal probabilities, and one sample is 
randomly taken from each interval. When considering implicit models, finite element analysis has been combined with 
probabilistic design methods in order to allow uncertainty in the system parameters to be taken into account. Several 
types of uncertainties can be distinguished: 

e Physical uncertainty: For example: material properties, dimensions, loading, etc.; 

e Statistical uncertainty: Due to limited sample sizes, probabilistic model is considered uncertain; 

e Model uncertainty: For example, mathematical models, numerical approaches due to simplifying 


assumptions, unknown boundary conditions, and unknown effects of other variables that are not included in the model. 


2.3.2 Reliability index 
The reliability is the inverse of failure probability. The relationship between the reliability R is related with 
the failure probability, it can be written as follows: 
P.=1-R (2) 
To estimate the reliability index, several techniques have been developed during the last five decades, namely 
FORM (First Order Reliability Methods), SORM (Second Order Reliability Method) and simulation techniques [11, 
14]. In FORM approximation, the probability failure is simply evaluated by 


P, x O(-B) (3) 


where (e) is the standard Gaussian cumulated function calculated through the following integral: 


O(Z) = ode. (4) 


l Z 
Ta! 
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For practical engineering applications, Eq. (3) gives sufficiently accurate estimation of the failure probability. 


In general, the nuclear and spatial studies necessitate very small values of failure probability, the failure probability 


should be: P, €[I 0° —10*] that corresponds to a reliability index B €[4.75—5.6] when using equations 3 and 4; while 


in structural engineering studies, the failure probability should be: P, €[1 0° —10°] that corresponds to a reliability 
index B €[3—4.25] (A detailed study on target safety indices can be found in Jeppsson [15]. In the next section, a 


simple 2D implant-bone study of solid and IAM stem designs 1s presented with account for statistical uncertainties. 

3. Results 

The Probabilistic Design System in ANSYS software analyzes a component or a system involving uncertain 
input parameters. The input parameters concerning geometry, material properties, boundary conditions, etc., are defined 
in ANSYS software. The variation of these input parameters is considered as random input variables and is 
characterized by their distribution type (normal, lognormal, etc.), as well as their distribution parameters (mean values, 
standard deviation, etc.). The important responses are defined as random output parameters. During a probabilistic 
study, the software executes multiple analysis loops to compute random output parameters as a function of a set of 
random input variables. The values for the input variables are generated randomly using Monte Carlo simulation. 

The studied model has six parameters that are regarded as random input parameters. Thus, six sources of 
E 


Cor ? 


uncertainty were considered in the present investigation: FE 


—_ Ev >Vcan>V cor ANd v,,. These sources are Young's 
modulus and Poisson's Ratio for the three studied layers (cortical, cancellous and metal layers). According to several 
references in the literature [16, 17], the maximum and the minimum values of Young's modulus are respectively 19.7 
GPa and 10.4 GPa for the cortical tissue, 5.6 GPa and 192 MPa for the cancellous tissue, and 120 GPa and 100 GPa for 
the titanium alloy. And the maximum and the minimum values of Poisson's ratio are respectively 0.33 GPa and 0.3 GPa 
for the cortical tissue and for the cancellous tissue, while they are: 0.36 and 0.3 for the titanium alloy. It is considered 
that the input parameters follow the uniform distribution law. To assess the accuracy of the results, it is performed with 
a high confidence interval of 99.87% (reliability level). The failure probability is noted to be 0.13% which corresponds 


to a reliability index equal to B =3 (Equations 3 and 4). Using Sampling Method, Monte Carlo based simulations were 


run, with 30000 simulations for the third loading case (L3) on the solid and IAM stems. A sensitivity analysis was 
performed to assess the influence of each parameter on the maximum von Mises stress values for the cortical and 


' ,o. and o” ). The evaluation of the sensitivities was based on both Rank- 


max ? max max 


cancellous bone, and for the metal (o 


order correlation coefficients between the input parameters and the output parameters. 


Table 2 
Statistics of the random input and output parameters for solid stem 


Standard a 
Parameter Mean ae Skewness Kurtosis Minimum Maximum 
Deviation 


E 


Can (MPa) 2896 1561 -1.73 x 107 6.65 x 10° 192.1 5600. 


E,, (MPa) 1.10x 10° 5774 2.36x 107 6.65 x 10° 1.00x 10° 1.20x 10° 
oe 0.3150 8.66x 10° -5.96x 10° 6.65x 10° 0.3000 0.3300 


Ec, (MPa) 1.51x 10° 2685 2.06x 107 6.65 x 10° 1.04 x 10° 1.97x 10° 


Tables 2 and 3 show the statistical results of the random input and output parameters for the solid stem and for 





the IAM stem, respectively. 
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Table 3 


Statistics of the random input and output parameters for the IAM stem 


Fig. 3 shows the probability density function histograms of output parameters. Fig. 3 a and b show the 





histograms of the maximum von Mises stress value (c.,,. ) for the solid stem, and the IAM stem, respectively. Fig. 3 c 
and d show the histograms of the maximum von Mises stress value (o~,.) for the solid stem, and the IAM stem, 


respectively. Fig. 3 e and f show the histograms of the maximum von Mises stress value (o””,.) for the solid stem, and 


the IAM stem, respectively. The probability density function calculates an appropriate number of classes based on the 
number of samples. The number of classes is equal to the number of bars which are presented in the histogram. The 
distance between the smallest and largest sample value is divided into classes of equal width. A histogram is modeled 
by counting the number of hits in the individual classes and dividing this number by the total number of samples. Thus, 
a histogram represents the relative frequencies of the random quantity that is drawn for. 
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Fig. 3. Histograms of the output parameters: 


(e) solid stem, and (f) IAM stem. 
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Fig. 4 shows the sensitivity evaluation of the output parameter with respect to the input random variables: 
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stem, and the IAM stem, respectively. Fig. 4 e and f show the sensitivities of the maximum von Mises stress value 


( 


2 
max 


stem, respectively. Fig. 4c and d show the sensitivities of the maximum von Mises stress value (o 


is 
o 


me ECO 
NUCOR 


significant parameters and insignificant 
i 


Significant: 
P ieee 
ECOR 
HUCOR 


) for the solid stem, and the IAM stem, respectively. The sensitivities of a certain random output parameter are 
a 
a 


one 
parameters. The sensitivity plots include only the significant random input parameters. 


modeled. The random input parameters are separated into two types: 


Insignificant: 


Insignificant: 


EM 
NUCAN 
NUM 


EM 
NUCAN 
NUM 


Significance level: 


2.500% 


Significance level: 


2.500% 








Output Parameter SMAX1 


‘ 
wo 
nM 
wu 
p 
; 
HM 
fr 
hy 
re) 
3 
oy 
rf) 
3 
oO 


> 


500% 


Significant 

a a 

a 
Insignificant: 
Significance level: 


EM 
NUCAN 
NHUCOR 
NUM 


° 
— 





Output Parameter SMAX2 


4) 
ol 


significant 
Significance Level: 


2.00 


Significant 
a 

Insianifi 
ECAH 

Ew 

HOCAH 
HECOR 

He 


ameter SMAx? 





sorUByoo| 


d) 


c) 


Par 


Output 


227i 


Advanced Engineering Research 2020. V. 20, no. 3. P. 216—224. ISSN 2687-1653 





1 i 


Significant: Significant: 
a) gg ™ -§ a 
a Insignificant: 2 Insignificant: 
; ECAN ECAN 
ECOR = ECOR 
aa EM _ EM 
NUCAN NUCAN 
-€ HUCOR 6 NUCOR 
5 Significance level: -5 — level: 
2.500% 2.500% 
4 
4 
oa 
oa 
a2 
ad 
= | 
ok 
a 


Output Parameter SMAXM 





Output Parameter SMAXM 


e) PD 
Fig. 4. Sensitivities of the output parameters: o| for (a) solid stem; (b) [AM stem; 07 ; (c) solid stem; (d) IAM stem; 


max max ? 


and oO” for (e) solid stem; (f) LAM stem 


max 


Table 4 shows the correlation coefficients between the input and output parameters for the solid and IAM 
stems. Here, it is shown the statistical interdependence between the input and output parameters. The values being 
closer to zero show that the two parameters are weakly correlated. However, the values being closer to | or -1, show 


that the two variables are highly correlated either in positive or negative sense, respectively. 
Table 4 
Correlation coefficients between the input and output parameters for solid and IAM stems 


Paramet 
ae (MPa) | of, (MPa) | o” (MPa) | o|,,(MPa) | 0%, (MPa) max (MPa) 


0 max 0 max 


or | oe | Ae | Ue ee ee 
Fen (MPa) 
ee ee ee ee ee oe 
ee eee ee eee ee 


4. Discussion. In this work, a probabilistic design strategy is established to compare in details the role of 
different input and output parameters considered when designing the solid and hollow (IAM) stems. Monte Carlo 
technique is used as a robust tool with a big number of simulations to provide with accurate results. The histograms of 
the probability density function of three output parameters are presented. The three output parameters have three 





different distribution forms: The maximum von Mises stress values for the cortical tissue (o' ) has the lognormal 


max 


distribution form. The maximum von Mises stress values for the cancellous tissue (o~, ) has the exponential 


max 


distribution form. The maximum von Mises stress values for the metal (o”) has the uniform distribution form. The 


skewness values of the input values (asymmetry) for the IAM stem differ from those for the solid stem, while there is 
no big difference when considering the output parameters. The sensitivity analysis for the output parameters with 
respect to the input random variables is next carried out in order to determine the input parameter influence. Three input 
) have different roles on the maximum von Mises stress values for the cortical tissue 


parameters (FE. ,E.. and v 


Can ? “~Cor 
1 
(o 


max 


Cor 


). Here, there is a small difference when comparing the diagram for the solid and IAM stems. For the maximum 


von Mises stress values for the cancellous tissue (o~,,. ), two input parameters (E,,. and E,,,,) play significant roles 


max 


when considering the IAM stem, while only one input parameter ( £,,..) when considering the solid stem. For the 


or 


maximum von Mises stress values for the metal (o™,. ), only one input parameter (-v,, ) has a significant influence for 


the solid and IAM stems. According to the correlation study, the correlation coefficients between the input parameters 
and the output ones for the IAM stem are much higher than those for the solid stem where several values are closer to 
zero. According to the presented probabilistic design strategy, the IAM stem has several advantages relative to the solid 
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5 Conclusion. A probabilistic design strategy is applied to find different probabilistic bounds with high 
reliability (confidence) levels. The results show that the IAM stem is much more advantageous than the solid stem 
especially. This study was limited to a 2D modeling in order to reduce the computing time since 30000 simulations 
were performed. However, in future work, it 1s recommended to perform a 3D modeling and to deal with anisotropy 
behavior for bone tissues. 
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